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The nonrelativistic annihilation of Majorana dark matter in the Sun to a pair of light fermions is 
chirality-suppressed. Annihilation to 3-body final states £^ f~V, where V = W, Z, 7, and £ and / are 
light fermions (that may be the same), becomes dominant since bremsstrahlung relaxes the chirality 
suppression. We evaluate the neutrino spectra at the source, including spin and helicity dependent 
effects, and assess the detectability of each significant bremsstrahlung channel at IceCube/DeepCore. 
We also show how to combine the sensitivities to the dark matter-nucleon scattering cross section 
in individual channels, since typically several channels contribute in models. 



I. INTRODUCTION 



Dark Matter (DM) particles X can become captured and trapped at the center of the Sun and the Earth. As the 
DM density grows over time, the accumulated DM can annihilate and produce a neutrino flux that is observable at 
a detector on Earth W. The annihilation channels and the final state decay products are determined by details of 
physics beyond the Standard Mode (SM). 

While it is often assumed that 2-body annihilations dominate, 3-body final states can be the leading contribution 
in models in which the DM candidate is a Majorana fermion. Such candidates arise, for example in supersymmetric 
models where the lightest supersymmetric particle is a neutralino. In such models, the cross section for dark matter 
annihilation to light fermions is severely suppressed [5]. 

The wavefunction of an initial state consisting of a pair of identical fermions must be totally antisymmetric, implying 
either L = 5 = 0orL = S' = l, where L and S are the orbital angular momentum and spin of the pair, respectively. 
For the latter case, the annihilation matrix element Ai{XX ff) is necessarily p-wave suppressed, and is thus 
proportional to w ^ 1, where v is the relative velocity of the DM particles. The s-wave initial state is CP-odd and 
has zero total angular momentum; if CP-violating effects are negligible, this state must annihilate to an L = 0, S' = 
final state. As the final state fermions // emerge back-to-back, they must possess the same helicity. Since particles 
and antiparticles of the same handedness arise from different Weyl spinors, s-wave XX — // annihilation requires 
that the initial state couple to both the /l and fa spinors, i.e., a mixed coupling to both L/R chirality. (For further 
elaboration of these issues see the Appendix of Ref. [3].) 

While fermion mass readily provides L-R mixing, it leads to a matrix element suppressed by mf/mx- If the 
mass term is the only source of helicity mixing, the 2-body XX — > // annihilation cross section is heavily suppressed. 
However, a 3-body final state containing an additional vector boson (VB) can be CP-even with vanishing total angular 
momentum, even if both fermions arise from the same Weyl spinor. As a result, the 3-body annihilation cross section 
is not suppressed by nij/rri^, and becomes significant despite the additional coupling factor (~ a). 

Radiative electroweak corrections to DM annihilation were recently considered for the gamma ray [IHS] , positron [3 
[8] and antiproton [8HIO] spectrum of the annihilations. Recently, solar DM signals from electroweak bremsstrahlung 
were investigated in Ref. In comparison, we consider each significant annihilation channel separately; the cor- 

responding event rates can be summed using annihilation branching ratios which depend on the details of a specific 
model. We also consider DM annihilation to left-handed and right-handed fermions separately. This is important, as 
the shape of the neutrino injection spectrum depends significantly on the helicity of the fermions (and in particular 
on their decay spectra). We also numerically propagate the neutrinos through the Sun and vacuum, with oscillations, 
scattering and r-regeneration fully simulated. 

In Section [llj we describe the model we adopt and compute the doubly differential 3-body annihilation cross 
sections. The injection spectra are presented in Section [lll[ and a description of neutrino detection in Section [IV| In 
Section |V) we investigate the discovery potential of the annihilation channels individually and in combination at the 



IceCube/DeepCore (IC/DC) detector. We conclude in Section VI 



2 



II. CROSS SECTIONS 



Here we briefly discuss the annihilation cross section in the case of 5?7(2)-singlet Majorana fermion dark matter 
X, with a Lagrangian similar to that of Ref. 12J, where X couples to SM fermions / through Yukawa terms, 

iint - VlXPlIvl + vIJPrXiiI + VRXPRfrtj,^ + yyPLXr,*j,^ , (1) 

where the yL,R are Yukawa couplings, r]L is a spin-0 SU{2) doublet and rjR is a spin-0 SU{2) singlet. Since X is a 
gauge-singlet under the SM, the two-body annihilation XX VV (where y is a vector boson) does not occur at 
tree level, even if kinematically allowed. In general, one may also write a mixing term ry^iyfl -f VrVl, whose coefficient 
is proportional to the Higgs vacuum expectation value. As L-R mixing lifts the suppression on XX — > // (see, for 
example, |13j). we restrict our attention to cases where such terms are negligible. 

Setting yR = without loss of generality, the leading contributions to the 2-body annihilation cross section are 
given by [5] 

,2 \ / „,4 .... 

■■''0{v^), (2) 

where r = rri^/m'^. As expected, the first (s-wave) term is suppressed by rnj/mj^ <C 1, while the second (p-wave) 
term is suppressed by w^. Note that at freeze-out v ~ 0.2 is not negligible, and the speed of DM particles in the solar 
core is much smaller than that in the galactic halo, v ^ 10^'^. 

In comparison, with the emission of a VB, the cross section for XX — !• ffV can be expanded as [Uj 



(3) 



As expected, the first two velocity-dependent (p-wave) terms are small compared to the corresponding term in the 
2-body cross-section because of the extra coupling factor g'^/in. However, the third (s-wave) term is velocity- 
independent and can be significant if r is not too large. In fact, for typical halo velocities, this 3-body annihilation 
process dominates the 2-body process for r < 0(10) 14j. Note that the s-wave term is not suppressed by the mass 
insertion as the VB spin cancels the total spin of the two fermions, which can have the same chirality. 

To compute the 3-body annihilation cross section, in addition to the interaction Lagrangian in Eq. (n]), we need 
the matter-gauge boson interaction vertices. These can be derived from the Lagrangian kinetic terms (m standard 
notation), 

, i (W+T^ + W-T ),,, ^3-sin^%Q , 



and 



,2 , 



|e ^ iW+T+ + W-T^) + 73 sin e^Q ^ ^^1 \ _ ^^2, ,2 

t V2sm^0w sm9wCos6w J J 

' -iW-T, + iy-T_) + e^J^-^^Z, + eAA' . (4) 



\/2sm^ 9w sin e'H' cos 6*1^ 

The matter-gauge boson interaction terms permit the tree-level process XX ffV to proceed via t- or w-channel 
exchange of t/l^r, with the VB radiated from either the external legs (i.e., final state radiation), or from the virtual 
ry (referred to as internal bremsstrahlung) . 
The total amplitude can be written as [2] 

iM = iKMa + M'r) + [Mb + Ms"") + {Mc + M§"^)] , (5) 

where the subscripts A-C refer to the diagrams in Fig.[T] We take the fermion masses to be negligible compared to the 
mass of the DM particle (nif jmx ^ 1). Using the FeynCalc 15] package, we find that the squared matrix element 
(summed over polarizations and averaged over initial spins) for annihilation to e.g., fhfhZ is 

_ (1 - 2 sin^ \Vl\'' (^(^i + :^2 - 1) + ^) (2 {xj - 2x, + - 2x^ + 2) - ^) 

\^\xx^S,hz- gsin^^vFCOS^^w m\{\ - 2x^ - rY{\ - 2x2 ^ r)^ ' ^ ' 
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FIG. 1. Feynman diagrams for the 3-body annihilation process, XX — > ffV. The incoming DM particles X have momenta ki 
and fc2. /, / are the outgoing fermion and antifermion with momenta pi and p2, respectively, rj is the mediator particle, and 
V is a gauge boson (7, Z, or W) with momentum k. 



where we use the notation of Ref. |5] and define xi — Ef-^/mx, 2^2 = Ef^/mx and 2:3 — Ez/mx, in the static center 
of mass frame with xi + X2 + X3 = 2. 

We separately compute the annihilation cross section to final states with any choice of fermion helicities. The 
differential cross sections, 



da 



"'■'^'d^Id^^ " 128^' ^ ' 

(where Vrei = vi — V2), for all final state channels and helicities are given in Appendix [Aj our results agree with those 
in Refs. [H [5J [101 E] . We have checked that integrating the differential cross sections in the mz — >■ limit, yields the 
results in Refs. [SI [71 |S1 [T3] . An analytic expression for the total cross section is given in Ref. [S]. 



III. NEUTRINO SPECTRA 



In this section we discuss neutrino injection from the leading annihilation channels with VB-bremsstrahlung. We 
focus on the couplings of DM to leptons because these are the most relevant to searches at neutrino detectors. In 
particular, lepton couplings produce neutrinos directly as part of the 3-body final state, and can provide a substantial 
contribution to the neutrino spectrum at high energy. However, lepton couplings do not contribute significantly to 
the DM capture rate; although DM can scatter off electrons in the Sun, such collisions do not result in DM capture 
because the momentum transfer is very small (since <^ mx)- Additional interactions between DM and Hght quarks 
thus provide the dominant contribution to the capture rate. It is worth noting that, although DM-quark interactions 
can also induce annihilation, these are unlikely to produce energetic neutrinos. These annihilation processes do not 
directly produce neutrinos in either the 2-body or 3-body final state, and the outgoing light quarks hadronize and 
stop before decaying, resulting in a very soft neutrino spectrum. Henceforth, we simply assume that there are some 
additional DM-quark interactions responsible for DM capture in the Sun. 

In general, the different annihilation channels are not independent of each other, and their branching fractions 
are determined by the couplings i/l.r (assuming the rj^ji states have degenerate mass). We present the spectra for 
individual channels, however, in order to illustrate which channels provide the hardest neutrino spectra. Moreover, 
for models in which the degeneracy of the 77L,fl is broken, the total injection spectrum can be found by summing the 
spectra of the individual channels after an appropriate rescaling. 

Lepton chirality plays an important role for both the annihilation cross section and the shape of neutrino spectrum 
arising from lepton decays. We investigate a scenario of flavor-independent lepton couplings, and a pure third- 
generation coupling (100% r) scenario. 
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FIG. 2. Normalized fermion spectra (solid) and VB spectra (dotted) for XX — )■ f fV for mx ~ 100 GeV and = 105 GeV 
from Eq. ([sjl. For light DM, off-shell VB emission is significant, resulting in 4-body spectra that deviate from these on-shell 
calculations. 

We calculate the tree-level matrix element for the 4-body annihilation, XX — > ff{V — > //) for all choices of 
fermion helicity. Thus, contributions in which the gauge boson is produced off-shell are included. The SM quantum 
numbers of rj^ ji are determined by gauge-invariance, and we assume that all scalar partners til^h share a universal 
mass niri = \jTmx- Consequently, the matrix elements are entirely determined by y^.R, f and mx- 

The leading contributors to the neutrino spectrum are 

• Primary neutrinos produced directly from the annihilation {XX — )• vvZ, l^DW^ , l^vW^). 

• Secondary neutrinos produced from the decay of primary l^v, Z ^ vv produced in the annihilation 
{XX ^ vDZ,l+l-Z,l-vW+,l+vW-). 

• Neutrinos from the decay of primary r/f 's produced in the annihilation process {XX — > fT{Z, 7), DrrW'^, fVrW~). 
We study the following channels: 

(1) XX ^ TLTLl 

(2) XX ^ TRTR-f 

(3) XX — > vv{Z — > vv) for z/e : t'p : i^r = 1 : 1 : 1 

(4) XX VrVr{Z VV) (100% VrVrZ) 

(5) XX ^TLfL{Z ^vv) 

(6) XX tbtr{Z ^ vv) 

(7) XX TLVr{W- l~Vl) + C.C 

(8) XX Ilv{W~' l~vi) + C.C for i/g : z/^ : i^T- = 1 : 1 : 1 

(9) XX l+l-{Z vv,fT), I = e,/x 

The subscripts L, R refer to the helicity of the Weyl spinor {e.g., is the antiparticle of a left-handed r~, which 
is a right-handed t+), and 'c.c' denotes the CP conjugate process. In addition to the three sources of neutrinos 
described above, for channel (9) we also include tertiary neutrinos arising from the decay of secondary r's produced 
from primary Z decay. This is the only channel for which the tertiary contribution is significant. 

Since the dense solar medium readily absorbs electrons and muons, channels (1-2) are the only channels with 
photon-bremsstrahlung that yield a significant neutrino flux. In these channels neutrinos arise from r decay, so the 
neutrino injection spectra are rather soft and are dominated by the flavor. Also, since the photon does not decay, 
the primary lepton and VB spectra can be obtained from the doubly-differential 3-body {XX — >■ ffV) annihilation 
cross section, 

da _ ^ da da _ ^ da 

dxi J^- ^ dx\dx2 ' dxy Jx- ^ dxidx2 

where xf = ^(2 — xy ± \/ Xy ~ ry) , x'2 — \ — x\— ry /A, = 1 — ry /{A{1 — xi)), and ry = {my j mxY' ■ The lepton 
and VB spectra are plotted in Fig. [2] for mx = 100 GeV and m,, = 105 GeV. 

Channels (3-6) involve Z-strahlung. Primary vv are produced in channels (3-4), and provide the dominant contribu- 
tion to the neutrino spectrum. Channels (5-6) lead to softer spectra because the primary fermions are r leptons. The 
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FIG. 3. Modification of tfie spectra from virtual VB production in tfie XX — > vuZ cliannel for mx = 50, 100 and 1000 GeV. 
The primary (solid) and secondary (dotted) neut rino spectra are sfiown. The 'On Shell' and 'w(ith) virtual' curves are generated 
using Calchep. Spectra obtained from Eqs. (All and ([8| are also shown. 



case of Z-strahlung in which the primary fermions are e/fi yields a still softer neutrino spectrum, and is treated sepa- 
rately in Channel (9). Channels (7-8) involve W-bremsstrahlung (with either 100% i^r or flavor-independent couplings), 
in which both primary and secondary neutrinos and charged leptons appear. The coupling in VK-bremsstrahlung is 
solely left-handed. We only consider leptonic W/Z decays, as the hadronic decays suffer from absorption in the solar 
medium and produce considerably softer neutrinos. 

To incorporate virtual W/Z contributions in Channels (3-9), we use the numerical package CalcHep [16| to compute 
the primary and secondary fermion spectra separately, and apply helicity-dependent decay to the t leptons, if present. 
In Fig. |3j we show the primary and secondary neutrino spectra from XX — )■ vDZ for mx = 50, 100, 1000 GeV, 
assuming either an on-shell VB, or allowing the VB to be off-shell. For DM mass above 100 GeV, the virtual 
contribution becomes subdominant and the primary spectra are reasonably well-described by 3-body cross sections. 

We display the neutrino injection spectrum from each channel for mx = 100 GeV, m,, = 105 GeV (Fig. |4]), and 
for mx — 1000 GeV, m,, — 1050 GeV (Fig. [s]). The antineutrino injection spectra are identical. 

Although the 3-body differential cross section for XX — TfifuV can be obtained from that of XX — )■ tlTlV by 
simply rescaling by coupling factors (see Appendix |A]) , the resulting neutrino spectra are quite different. This can be 
seen by comparing the spectra for Channels (1) and (2), and for Channels (5) and (6) in Figs. [4] and [5) The reason 
for this difference is that the neutrino spectrum arising from the decay of a r depends on its helicity; the neutrino 
spectrum arising from a highly-boosted left-handed r differs markedly from that of a right-handed r. For each channel, 
the fraction of left/right handed t is set by the couplings. W decay produces 100% left-handed r while the left /right 
fraction from Z decay is 57%/43%. See Ref. [T7] for a description of our treatment. Note that the channels dominated 
by r decay lead to a rather soft neutrino spectrum; its power-law shape is less distinctive compared to channels that 
yield primary neutrinos. 

Due to the finite mass of the W/Z bosons, the neutrino spectra depend non-trivially on mx- Consider, for example, 
the XX — > lyDZ channel: if mx ^ mz/2, then the cross section for producing an on-shell Z is suppressed by the 
phase space of the 'primary' neutrinos (Ei, ~ mx — mz/2). The primary neutrino can be much softer than mx, even 
below that of the secondary neutrinos, as shown in Fig. |3] As can be seen from Fig. |5j for mx ^ mz, the vDZ and 
liyW channels are by far the hardest. 

It is worth noting, however, that the finite mass of the the W and Z gauge bosons can also enhance the annihilation 
cross section to those channels. The reason for this enhancement is that for some regions of phase space, the final 
state can only have vanishing total angular momentum if the gauge boson is helicity-0. These contributions must 
vanish in the my limit, implying that these terms in the squared matrix element scale as my/m^. 

The mediator mass also has a noticeable impact on the shape of neutrino spectrum. In general a larger 
leads to a neutrino energy distribution that is less peaked at the end-point; this effect could slightly enhance the signal 
rate, since at lower energy (yet above detector thresholds) neutrinos suffer less attenuation from scattering. However, 
for the values of mx under consideration, varying m,, yields only insignificant changes to the shape of the spectrum. 
Moreover, the annihilation cross section is suppressed by l/r^. We use a low mediator mass r = 1.1 throughout. 
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FIG. 4. i/^ and i/t injection spectra from each annihilation channel (as labeled), including r decay contributions, for mx = 
100 GeV and = 105 GeV. The Ue injection spectra are identical to that of Vfj,. For comparison, the gray dotted curves show 
the primary lepton spectrum in each case. The antineutrino injection spectra are identical. 



IV. NEUTRINO DETECTION 



After injection at the solar center, the neutrinos propagate through the Sun, and then through vacuum to Earth. 
In doing so, they oscillate and undergo scattering. We take the neutrino mixing parameters to be 

Smj^ = 8.1 X 10"^ eV^ , Sml^ = 2.2 x 10"^ eV^ , 6*12 = 33.2° , 023 = 45° , 6*13 = 8.8° , 

where drn'^^ = — and we use the value of ^13 recently measured by the Daya Bay experiment jLSj. For details of 
our simulation of neutrino propagation (including oscillations, tau-regeneration, and energy losses due to collisions), 
and muon detection at IceCube/Deepcore, see Refs. pTl [TO] . 

To assess the IC/DC sensitivity we choose a muon energy window, Eth < < mx, where Eth is the detector 
threshold. We target a 3ct detection, TV'^ — 3\/iV°*™, where A^^^f and iV"*™ are the number of signal and atmospheric 
background events, 

N^^3= dt V{6{t),E^)---dE^, (9) 
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FIG. 5. Similar to Fig.|4) for mx = 1000 GeV and = 1050 GeV. 



dt 



V{e{t),E,) dE,. 



Et 



(10) 



Here, and (J)"*™ are the muon fluxes generated by the signal and atmospheric neutrinos incident on ice, respec- 
tively, and V is the effective dimension of the IC /DC detector relevant to the event type. For neutrinos arriving from 
the direction of the Sun, most charged leptons scatter within a cone of half-angle. 



59 = 20° 



10 GeV 

mx 



(11) 



see, e.g., |20J . Equation (11) gives the intrinsic scattering angle for incoming neutrinos with energy ~ mx, which is 
comparable to the IC detector's angular resolution, and which we use to estimate the atmospheric background. Note 
that since iV°*™ scales quadratically with So, the constraint on the DM annihilation rate scales linearly with 5g. 

Besides the variation with zenith angle and event energy, V depends on whether the muon event is up-going 
or contained |21j . Up-going (contained) events refer to upward going muon tracks that start outside (inside) the 
instrumented volume of the detector. A discussion of the effective detector dimensions for these two types of events 
can be found in Refs. [171 121| . In this analysis we consider up-going muon events for IC and contained events for 
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FIG. 6. The Sa sensitivity to aff and aff of IC up-going, IC contained, and DC contained events, with one year of data. 
We assume that the DM capture and annihilation processes are in equilibrium, and set the branching fraction to each channel 
equal to unity. 



DC, and a 1 km^ effective volume with a 70 GeV energy threshold as an estimate for IC contained events. The 
effective IC area for up-going events falls rapidly below Eth = 60 GeV and we assume Etu — 10 GeV for the DC 
subdetector. Since observations track the trajectory of the Sun, we use the zenith-angle dependent atmospheric 
flux measured by Super-Kamiokande \21^ . The DM-induced muon flux has an energy spectrum that is determined by 
the DM capture rate and the neutrino spectrum. The time integral in Eqs. (9p0| spans half a year for IC and a full 
year for DC, which has 47r angular coverage. Descriptions of our calculations of up-going and contained event rates 
are provided in Refs. [T71 119) . 



V. DISCOVERY POTENTIAL 

IceCube's Scr sensitivity to DM annihilation can be determined by setting N'^^^ — 3\/N^^. We assume that the 
DM capture and annihilation processes in the Sun are in equilibriumj^ in which case a constraint on the annihilation 
rate can be directly translated into a constraint on ctat x BFj , the product of the DM-nucleon scattering cross section 
and the branching fraction to the annihilation channel in question. We normalize the branching fraction to leptonic 
channels i?F(leptons) to unity, so that J^j = i3F(leptons) = 1. Therefore, the sensitivities to the DM-nucleon 
cross section presented below should be divided by the actual value of i3i^(leptons). In Table |l] we present the 
sensitivity to each annihilation channel (using each of the three event samples) for mx = 50, 100 GeV and 1 TeV. The 
3(T sensitivity to the cross section for spin-independent (SI) and spin-dependent (SD) scattering are listed separately as 



^ Whether dark matter capture and annihilation are in equilibrium is very model dependent. We pick the equilibrium case as a reference 
signal rate and present the rest of the calculation in as model-independent a manner as possible. Given the details of a specific model, 
the total annihilation rate (as a fraction of the capture rate) and the annihilation branching fractions can be determined, and our results 
can be translated into a bound on the particular model. 

Note that for r = 1.1, yL = Vr = 1 and mx = 100 GeV, the total annihilation cross section to leptonic final states is 2.4 pb, and 
for yi^ = yji = VStt and mx = 1000 GeV it is 0.8 pb. For annihilation cross sections of this order, and for scattering cross sections of 
the size to which IC/DC is sensitive, equilibrium will hold. There can also be a contribution to the total annihilation cross section from 
quark final states (which provide only a subleading contribution to the neutrino spectrum). 
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TABLE I. The 3a sensitivity to the spin-dependent (SD) and spin-independent (SI) DM-nucleon scattering cross sections for 
each annihilation channel with one year of data. The number of atmospheric background events jV"*™ in one year at DeepCore 
and 180 days at IceCube is also provided. The muon energy window runs from the experimental threshold (10 GeV for DC, 
60 GeV and 70 GeV for IC up-going and contained) to mx- 



af/ and crff , respectively. The dependence of the equilibrium annihilation rate on cr^v is calculated using the method 
described in Ref. '5(T. The corresponding background rates are also listed. In Fig. 6] we plot the 3cr sensitivity to the 
DM-nucleon scattering cross section (setting the branching fraction to each channel equal to 1), under the assumption 
that equilibrium holds between capture and annihilation in each individual channel. The combined sensitivity of the 



for each channel. 



three event samples obtained from ^ X]i=i J — 3, is shown in Fig. 7 

For any particular model, there will be several annihilation channels "and the total event rate will be the sum of 
the contributions from each channel. To compute the sensitivity to the DM-nucleon cross section from a combination 
of channels, cr^'"'', we utilize the fact that signal significance scales linearly with the cross section. We denote by 
cr^ j the 3cr sensitivity of IC/DC in annihilation Channel (j), using event sample i (either DC contained events, IC 
up-going events, or IC contained events), as reported in Table [ij Denoting the branching fraction to Channel (j) by 
BFj^ we find the 3ct sensitivity for event sample i from 



3 



(12) 



The 3a combined sensitivity from a combination of channels and all three event types is then 

1 

3 



comb 



.i=l ^N,t 



(13) 



In Table pi we show the 3a sensitivity to afP and af/ for mx = 100, 1000 GcV for three models: DM with couplings 
only to left-handed third generation leptons {yfj = 0), DM with couplings only to to tjj (j/^ = 0), and DM with equal 
couplings to left-handed and right-handed third generation leptons {i/l = Vr)- 
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model 


(7^° {mx = 100 GeV) 


af/" (mx = 1 TeV) 


af/ {mx = 100 GeV) 


af/ {mx = 1 TeV) 


yR = 

yL = o 

Vl = Vr 


1.4 X 10"^ pb 
3.4 X 10"^ pb 
1.4 X 10"^ pb 


1.3 X 10"^ pb 
2.1 X 10"^ pb 

1.4 X 10-^ pb 


5.3 X 10"* pb 

1.4 X 10"^ pb 
5.6 X 10"* pb 


1.2 X 10"** pb 
1.9 X 10"** pb 

1.3 X 10"** pb 



TABLE II. The 3a sensitivity to af/^ and af/ (with 1 year of data) for three models with Yukawa couplings only to third 
generation leptons: DM coupling only to left-handed leptons {yn = 0), DM coupling only to tr {i/l = 0), and DM with equal 
couplings to left- and right-handed leptons {yL = Vr)- 



VI. SUMMARY 



If dark matter is a Majorana fermion, its annihilation in the Sun to Standard Model fermions is both chirality and 
velocity-suppressed. Then, the annihilation may primarily be through 3-body processes, with the emission of a gauge 
boson. The neutrino spectra from such channels can differ dramatically from the spectra from 2 — 2 processes. 

Dark matter couplings to left-handed leptons necessarily open-up 3-body annihilation channels in which neutrinos 
are produced directly. Moreover, the branching fractions to these channels are usually large. The neutrino injection 
spectra are typically hard, providing for interesting detection possibilities at neutrino detectors. 

We considered a model in which S'f/(2)-singlet DM couples to SM leptons (either left-handed or right-handed) via 
exchange of a new scalar rj^^n. For this model, we calculated the 3-body differential annihilation cross sections and 
the neutrino injection spectra with a full treatment of the helicity correlations of the gauge boson and t lepton decays. 
We determined the muon event rates at IceCube/DeepCore arising from each annihilation channel, accounting for 
neutrino propagation effects, including oscillation, scattering and regeneration. 

We calculated 3a sensitivities of IC/DC to the DM-nucleon scattering cross section for several 3-body channels. 
The different channels are of varying utility in constraining dark matter models; channels with primary neutrinos lead 
to the best sensitivity. 

We then showed how to combine the sensitivities in individual channels, to obtain the sensitivity for a combination 
of channels as may arise in models. 
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Appendix A: Vector boson emission cross sections 

As shown in Ref . [6] , the analytic form of the differential cross section is independent of the emitted electroweak 
gauge boson. Consequently, the set of differential and total cross sections for the various channels can be obtained 
from the following equations: 



da 



Vrel 



dxidx2: 



or3 



e2 (1 - 2 sin^ ewf IvlI^ (^(xi + - 1) + ^) (2 {xj ~ 2x, + xj - 2x2 + 
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XX^fLfLZ 



10247r3 sin^ dw cos^ 9w 



m\{l - 2xi - r)2(l - 2x2 - r) 



where [S] 



(7 



4 sin dw 

. 2 „ \2 ^XX^fL h ZIVL^VR 



(1 - 2sin^ dw) 
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